Localized states and interaction induced delocalization in Bose 
gases with quenched disorder 



o 
o 

oo 



I 

O 

o 



(N 
> 

in 

00 

o 

00 

o 



X 



G.M. Falco^, T. Nattermann^ and V.L. Pokrovsky^-^ 



Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Str. 77 D- 50937 Koln, Germany 

Dept. of Physics, Texas A&M University, College Station, TX 77843-4242 

Landau Institute for Theoretical Physics, Chernogolovka, Moscow District 142432, Russia 



PACS 03 . 75 . Hh - 03.75.Kk 

Abstract. - Very diluted Bose gas placed into a disordered environment falls into a fragmented 
localized state. At some critical density the repulsion between particles overcomes the disorder. 
The gas transits into a coherent superfluid state. In this article the geometrical and energetic 
characteristics of the localized state at zero temperature and the critical density at which the 
quantum phase transition from the localized to the superfluid state proceeds are found. 



■ Introduction. — The interplay between interaction 

■ and disorder is an important paradigm of condensed mat- 
' ter physics. In 1958 Anderson [1] showed that in disor- 
. dered solids a non-interacting electron may become lo- 

■ calizcd due to the quantum interference. A phenomeno- 
logical theory of locahzation [2, 3] concluded that non- 
interacting electrons in one and two dimensions are always 
localized. In three dimensions the localized and extended 

^ states are separated by the mobility edge. States with 
energy significantly below this edge in 3 dimensions are 
strongly localized. They appear in rare fiuctuations of the 
quenched random potential [4-6]. These instanton-type 
states broaden and eventually overlap with growing en- 
ergy. A system of non-interacting fermions in the random 
potential transits from the insulator to metal state when 
its Fermi energy exceeds the mobility edge. Thus, the 
Pauli principle delocalizes fermions in 3 dimensions, but 
leave them localized in lower dimensions. The common be- 
lief is that the repulsive interaction suppresses the localiza- 
tion. So far this problem was studied only in the limit of a 
weak disorder [7,8]. Therefore, the interaction induced de- 
localization transition remains beyond the frameworks of 
the theory. The metal-insulator transition in 2 dimensions 
was observed in experiments [9] suggesting the decisive 
role of interaction. The physical picture changes drasti- 
cally for bosons. The non-interacting bosons condense at 
a single-particle state with the lowest energy. In a homoge- 
neous system it leads to a coherent quantum state known 
as the Bose-Einstein condensate (BEG). Examples are su- 
perfluid phases of He [10], superconductors [11], BEG of 



ultra-cold alkali atoms [12, 13] and of excitons in semicon- 
ductors [14]. BEG still persists when a small amount of 
disorder is added to the system. BEG in a random en- 
vironment was observed in the superfluid phase of '*He 
in Vycor glass or aerogels [15] , in "^He in aerogels [16] 
and in ultra-cold alkali atoms in disordered traps [17-23]. 
But in a random environment and in the absence of in- 
teraction, all Bose-particles fall into the lowest localized 
single-particle state. Such a ground state is non-ergodic 
since its energy and spatial extension depend on a specific 
realization of the disorder. An arbitrary small repulsive 
interaction redistributes the bosons over multiple poten- 
tial wells and restores ergodicity. Hence, contrary to the 
fcrmionic case, the perturbation theory with respect to 
the interaction strength is invalid. At low temperature, 
the Bose systems display superfluidity provided the den- 
sity n of bosons exceeds a critical value ric- At either 
weak disorder or strong interaction, i.e. at n ^ n^, the 
disorder corrections to the superfluid density rig (and the 
condensate density no) are small [24-26] . These correc- 
tion blow up with the interaction decreasing, signaling the 
breakdown of the theory. 

We present an alternative approach to the problem 
of the interaction-induced delocalization starting from 
deeply localized state of the Bose-gas in a random poten- 
tial. We present a simple and visual picture of the deeply 
localized state, which decays into remote weakly coupled 
fragments. We give a geometrical description of fragments 
and their distribution in space. At a critical density Uc, 
which we express in terms of the disorder characteristic 
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and interaction strength, the increasing tunneUng of par- 
ticles between fragments leads to transition from the ran- 
dom singlet state to the coherent superfluid. 

Single-pcirticle levels in an uncorrelated random 
potential. — The random environment produces a ran- 
dom potential ?7(x) for the bosons. We assume that ?7(x) 
is Gaussian distributed with zero average and short range 
correlations 



([/ (x) [/ (x')) = k'(5(x - x') 



(1) 



We will consider briefly the long range correlated case at 
the end of this article. In the absence of interaction the 
single-particle wave functions obey the the Schrodinger 

equation 

|-VV + (S-t/(x))^-0. (2) 

Its energy levels E \U (x)] are functional of the potential 
U (x) . The only characteristic of the random potential 

K together with the Planck's constant Ti and the mass m 
establishes the scales of length and energy: 



(3) 



which we call the Larkin length and Larkin energy, re- 
spectively [27]. The density of states v{E) belonging to 
(2) in the limit E < Q,\E\ :^ E was calculated in [4-6] 
(for a complete summary see [28]). For 3d system with 
the volume Q. it reads: 

v[E) = ^{6{E-E[U (x)])) ~ 7\^(£)e-(l^l/^)'^', (4) 

where we absorbed a numerical constant in the exponent of 
viE) in the definition of £. As we show below the precise 
form of the prefactor Af{E) is not relevant for our consider- 
ation. In a large 3d volume the states with energy E ^ £ 
are delocalized, whereas the states with negative energy 
sufficiently large by modulus E < Q and \E\ ^ £ are 
strongly localized. The threshold of localization is a posi- 
tive energy of the order of £ [29] . In the interval between 
£ and —£ the transition from the extended to strongly lo- 
calized states proceeds. The latter are supported by rare 
fluctuations of the random potential, which form a poten- 
tial well sufficiently deep to have the negative energy E as 
its only bound state. Let us introduce the spatial density 
nyj{E) with the energy less than E. It is related to the 
DOS by equation nw{E) = j^^v{E)dE. For deep levels 
it can be also considered as the spatial density of states 
nu,{R) with the radius less than R, where R = li/ y/2m\E\. 
For such states (R) is proportional to a small exponent 
exp(— -^/j E \/£) = exp(— £/i?). From the dimensionality 
consideration it follows: 



m = R~'f 



exp - 



R 



(5) 



The function f{x) can be found from Ref. [30] to be pro- 
portional to f{x) ~ cc" with a = I. It will be inessential 



for further calculation. The average distance d{R) be- 
tween the wells of the radius less than R reads: d{R) = 
riw^^^ = i?/~^/^exp (^). Thus, the distances between 
the wells are significantly larger than their sizes. The tun- 
neling factor t (R) between two typical wells with the ra- 
dius R of the same order of magnitude is given by a semi- 
classical expression t {R) — exp (— ^ / \p\ dl), where the 
path of integration connects the two wells. By the order 
of magnitude p ^ h/R and the length of the integration 
path is r^d{R). Thus, ^ / \p\ dl « d/R w /-^/^exp (^) 
and 



t {R) = exp 



-1/3 



exp 



3RI 



(6) 



At i? ~ £/3 or ~ —9£, the distances between the opti- 
mal potential wells become of the same order of magnitude 
as their size R. Simultaneously the tunneling amplitude 
between the wells becomes of the order of 1. The poten- 
tial wells percolate and tunneling is not small, but the 
states still are not propagating due to the Anderson local- 
ization [1]. 

Bose gas in a large box with an uncorrelated 
random potential. — In the ground state of an ideal 

Bose gas in a large box with the Gaussian random po- 
tential all particles are located at the deepest fluctuation 
level. In the box of cubic shape with the side L the deep- 
est level which occurs with probability of the order of 1 
has the radius R determined by equation: L^riu, (R) = 1, 
i.e. R ^ 3in(L/£) • ^'^^ prefactor / introduces a negligible 
correction to the denominator of the order of In(ln^). 

The corresponding energy is ~ — 9f (in^)^. As we 
already mentioned such a state is non-ergodic since the 
location and the depth of the deepest level strongly de- 
pends on a specific realization of the disordered potential. 
Therefore, the average energy per particle and other prop- 
erties averaged over the ensemble has nothing in common 
with the properties of a specific sample. Even an infinitely 
small repulsion makes the system ergodic in the thermo- 
dynamic limit, i.e. when first the size of the system grows 
to infinity and then the interaction goes to zero. In a 
sufficiently large volume any physical value per particle 
coincides with its average over the ensemble. The reason 
of such a sharp change is that, at any small but finite in- 
teraction, the energy of particles repulsion overcomes their 
attraction to the potential well when the number of par- 
ticles increases. They will be redistributed over multiple 
wells. Since the distribution of wells in different parts of 
sufficiently large volume passes all possible random config- 
urations with proper ensemble probabilities, the ergodicity 
is established. Below we find how the interacting particles 
eventually fill localized states. In a real experiment the 
Bose gas may be quenched in a metastable state depending 
on the cooling rate and other non-thermodynamic factors. 
This is what M.P.A. Fisher et al. [31] call the Bose glass. 
Such a state is also possible in the case of weakly repul- 
sive Bose gas. However, as it will be demonstrated later, 
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in the case of cooled alkali atoms the tunneling amplitude 
still remains large enough to ensure the relaxation to the 
equilibrium state in 10~^-;-10~^s. Our further estimates 
relate to the real ground state. As in the Bogolyubov's 
theory [32] we assume that the gas criterion na^ ^ 1 is 
satisfied. Here n = JV/fi is the average particle density; N 
is their total number and a is the scattering length. Im- 
plicitly our considerations takes in account the change of 
the optimal potential well due to the interaction. 

Let the Bose gas with the average density of parti- 
cles n fill all potential wells with the radii less than R 
in the ground state. The average number of particles 
per well is N^j (R) = n/riw (R)- The local density inside 
the well of the linear size R is np{R) — ^^^rP^- The 
gain of energy per particle due to random potential is 
E (i?) = — 2mR^ ' ^^'^ repulsion energy due to interaction 
3h Af„(fl)a ^]^gj.g used the well- 



2ml 

is equal to gUp (R) 
known relation for an effective potential field induced by 
a gas of scatterers [33]. Minimizing the total energy per 

particle Etot (R) = + mfl^^^^° ^^^"^ ^ ^'^'^ 

the value of R corresponding to the minimum of energy at 
fixed n with the logarithmic precision: 



R{n) 



C 



ln{nc/n) 



(7) 



where 



= 3£2 



denotes the critical density. The 



factor / in equation (5) leads to corrections of the type 
ln(ln(nc/n)) which can be neglected. Further we put / — 
1. The distances between the filled wells according to the 
corresponding expression d{R) for single-particle states 
reads d{n) = £(ln(nc/n))~^ {uc/n)^^^. They strongly ex- 
ceed the average size of the potential well (7) at n ^ n^- 
At the same condition the chemical potential of atoms can 
be estimated as fi{n) — — ^L-n = — i fln^)^. The 



2mR^(r. 



tunneling amplitude t (n) between two wells separated by 
a typical distance d (n) can be found by employing the 
single particle result (6): 



t (n) — exp — [nc/n)^^^ 



(8) 



Thus, the Bose gas at n <^ ric is fragmented into multi- 
ple clusters of small size R (n) separated by much larger 

distances d (n) and containing about C/ 3a (in ^)'^ par- 
ticles each. The amplitude of tunneling between the wells 
depends on the scattering length in a non-analytic way and 
is exponentially small for weak interaction. Therefore, the 
number of particles in each cluster is well defined. As a 
consequence, the phase is completely uncertain. Such a 
state is a singlet with non-uniformly distributed particles, 
a random singlet: the ground state is non-degenerate. The 
compressibility = J In ^ is finite as expected for the 
Bose glass phase [31]. 

Bosons in atomic traps. — Our results can be easily 
extended to bosons in harmonic traps characterized by a 
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Fig. 1: Regime diagram of atoms in traps: uncorrelated disor- 
der. R denotes the size of the single existing atomic cloud. L 
is the size of the cloud of fragments. 



potential 



^trap — 



2m 



(9) 



where we introduced the oscillator length i — ^Jh/ {rauj). 
This section partly overlaps sligthly with our previous 
work [34]. The energy of the bosons includes now four 
terms: the kinetic energy, the confining potential energy 
of the trap, the repulsion from other particles and the 
energy of the random potential. Two of them, the inter- 
action with the trap and the random potential tend to 
confine and localize the particle. Going through essen- 
tially the same steps as before, we can distinguish four 
different regimes (Figure 1). 

1. Weak disorder and weak interaction: 3Na <ti i ^ C . 
In this case the interaction can be neglected. Minimizing 
the remaining terms, the kinetic energy and energy of the 
trap, we find R — £. Physically it means that all particles 
are condensed in the oscillator ground state. 

2. Weak disorder and strong interaction: £ ^ C, £ <^ 
3Na . Neglecting the kinetic energy and minimizing re- 
maining energy of traps plus the repulsion energy, one 
finds the result known as Thomas-Fermi approximation 



R 



[13]: 

3. Strong disorder and weak interaction : 3Na <^ C £ . 
In this range of variables the non-ergodic phase is realized. 
Since interaction is negligible, the particles find a random 
potential well with the deepest level and fall into it. Let 
such a well can be found at a distance ~ L from the trap 
center. Its depth typically is about 95 In^ (!//£). This 
gain of energy must be not less than the loss of the trap 
energy muP'LF' /2. A typical value of L appears when both 
this energies have the same order of magnitude. Thus, 
L « Q^/2{£'^/C)\n{£/C). A typical size of the weU is 
Rk CI (6 In {£/C)). 
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4. Strong disorder and moderate interaction: C <^ 
3Na <C £■ In this case the ergodicity is restored. Our 
experience with the gas in a box prompts that the gas 
cloud is spht into fragments each occupying a random 
potential well from very small size till same size R de- 
pending on N. The typical disorder energy per particle is 
H = —£ (in ^) . It becomes equal to the trap energy at 
the distance L ^ {l"^ / InF where F is a new dimension- 
less parameter 

1)6 ^ 

-. (10) 



HNa/i) , L^Cfm- ' ^ 



Therefore, the average density is n ^ NC^ /£^. The state 
of the Bose gas is fragmented and strongly localized when 
F is large; the transition to delocalized superfluid state 
proceeds when this ratio becomes ~1. The phase diagram 
is shown in Fig. 1. Note the counter- intuitive dependence 
of the size on the number of particles: the cloud slightly 
contracts with increasing number of particles. It happens 
because the number of particles in each fragment increases 
more rapidly with the average density than the number of 
fragments. 

Correlated disorder. — So far we considered un- 
correlated disorder (1). Our results can be extended to 
random potentials with a finite correlation length b and 
strength Uq = ^ (f/^ (x))- We quote here the results 
without derivation, which can be found in Ref. [36]. As 
long as 6 ^ £ = ^^^^-1^3 the results of the previous 
considerations remain correct. In the opposite case the 
optimal potential wells have the width h and, contrary 
to the short range correlated case, they contain many 
bound states (of the order of {bj Cf'^'^ ^1). It is con- 
venient to introduce a new length scale B — b{C/hf'^^ = 

(3/47r)^/''(;iVmC/o) ^ In the following we restrict our 
consideration to the case b ^ £, i.e. b ^ B. The den- 
sity of states in this case is i^{E) ^ exp[— {E / Uq)'^] [28]. 
The critical density is Uc ^ l/{aB'^) (this result has 
been found before in [35]) and the typical size of a frag- 
ment is i? w b{ln{nc/n))~^/^. The distance between frag- 
ments is d(n) w b{nc/ny^^ and the tunneling coefficient 
is t{n) « exp[-(&/-B)K/n)i/3]. 

In the case of a harmonic trap we again find four 
different regimes (Fig. 2). The relevant parameter is 
F = e^/{NaB^) w Hc/n. For F w 1 the transition to the 
superfluid phase proceeds. All results can be extended to 
lower dimensions [36]. 

Conclusions. — Four parameters can be controUably 
and independently varied in the experiment. They are: 
number of particles N; the frequency lu or equivalently 
the strength of the trap; the scattering length a (it can be 
varied by approaching one of the Feshbach resonances); 
the strength of disorder Uq. Using this freedom it is feasi- 
ble to pass all regimes described above. A simple estimate 
shows that, at 6 l^m, the transition from uncorrelated 
to strongly correlated regime proceeds at frequency of dis- 
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Fig. 2: Regime diagram of atoms in traps: correlated disorder. 
R denotes the size of the single atomic cloud, L denotes the 
size of the fragmented state. 



order potential ujd = y''2UQ/mb'^ ^ Ik Hz which is acces- 
sible. 

Simplest experiments are the measurements of the cloud 

size L as a function of different variable parameters in the 

regime of multiple localized fragments. Theory predicts 

that in the regime of uncorrelated disorder the size of the 

cloud is proportional to Uq /uj . It also predicts very weak 

dependence of the size on the number of particles ~ In N. 

In the case of strongly correlated disorder the size of the 

1/2 

cloud is proportional to luUq' ; the dependence on N also 
is weaker than in the uncorrelated regime: L cx (In A^)^^"*. 

It would be important to observe the transition from 
non-ergodic state with one or few fragments to the er- 
godic state with many fragments and check that it hap- 
pens at N = C/3a for uncorrelated disorder and at 
N = (b^/3aB'^^ for strongly correlated disorder. 

Another feasible experiment is the time-of-flight spec- 
troscopy after switching off both the trap and the random 
potential. In this experiment the distribution of particles 
over momenta (velocities) is measured. Its width Ap is 
associated with the average size of the fragment R by the 
uncertainty relation Ap = h/R. It gives the opportunity 
to check the equation R = £/lnF for the uncorrelated 
disorder or i? = £/lnF for correlated disorder. Installing 
a counter close to the trap, at a distance comparable to 
the size of the trap, would allow to register the oscillations 
of the particle flux due to discrete character of the frag- 
mented state. This is an opportunity to find the distances 
between fragments and compare theory with experiment. 

The transition between localized and delocalized coher- 
ent state in the random potential was found in several ex- 
periments (see Introduction). We propose to make more 
detailed measurement of the transition manifold and check 
our predictions. 
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An important question is whether the relaxation to the 
ground state can be roachcci during a reasonable time in- 
terval compatible with the time of experiment. We analyze 
this question for the uncorrelated or weakly correlated dis- 
order. In this case the relaxation time due to tunneling can 
be estimated as r = 2'KUj~^t~^, where u;„ ~ ^ (InT)^ is the 
characteristic frequency of the optimal potential well and 
t ^ cxp [— r^/^] is the tunneling coefficient (sec cq. (8)). 
For numerical estimates we accept T ~ 125, i ~ 10/xm, 
b'^ C-^ 1/xm, a ~ 0.01/xm, N ~ 27, 000. Then t''^ = 148 
and T ^ 0.06s. The Larkin length can be increased by 
decreasing the amplitude of the random potential. Simul- 
taneously, at fixed values A'', £ and o, the value F decreases 
as . This example shows that the equilibrium is ac- 
cessible, though it is difficult to reach large ratio £/b. 

The closest to ours was the approach developed in the 
work Lugan ct al. [20]. Apart from the fact that these 
authors considered only the 1-dimensional case, the main 
difference between our and their problems is that they con- 
sidered the random potential with the exact lower bound- 
ary [/b and with on-site distribution function W[U] oc 
exp[— c(f/— C/b)] instead of a Gaussian distribution. Such a 
distribution allows deeply localized states only at energies 
E close to the exact lower boundary Ui,. The correspond- 
ing fluctuations have the width R the broader the closer 
is E to Uh- It is clear that these levels arc very different 
from those discussed above. If the random potential has 
the exact lower boundary, our theory is valid only if this 
boundary is separated from the most probable value of 
the potential by an energy interval strongly exceeding the 
energy dispersion. Then the localized states of our theory 
appear at intermediate energies between dispersion and 
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